Tree-body loss of of trapped ultracold 87 Rb atoms due to a Feshbach resonance 
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The loss of ultracold trapped atoms in the vicinity of a Feshbach resonance is treated as a two- 
stage reaction, using the Breit-Wigner theory. The first stage is the formation of a resonant diatomic 
molecule, and the second one is its deactivation by inelastic collisions with other atoms. This model 
is applied to the analysis of recent experiments on 87 Rb, leading to an estimated value of 7 x 10 -11 
cm 3 /s for the deactivation rate coefficient. 
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The phenomenon of Feshbach resonance has received 
recently an increased attention due to its application to 
Bose-Einstein condensation (BEC) (see Ref. and ref- 
erences therein). Rs most outstanding effect is a drastic 
change of the elastic scattering length as the collision en- 
ergy of an atomic pair approaches the energy of a bound 
level belonging to another electronic or hyperfine state. 
The resonance can be tuned by applying an external mag- 
netic field, as has been proposed in Ref. @ in order to 
control the BEC properties. Applications include a con- 
trolled BEC collapse and bright solitons in BEC J 
as well as a formation of molecular BEC EIJH I 
an atom-molecule coherent superposition UflOl, and an 
entangled atomic gas 0. 

Another effect of the resonance is the abrupt increase 
in atom loss due to inelastic collisions of the resonant 
molecules 0, |M and to the formation of non- 
condensed atoms [21 HI IS- The determination of the loss 
parameters is important for an appreciation of the out- 
come of applications of Feshbach resonances. We present 
here an estimate of the rate coefficient for the deactiva- 
tion of vibrationally excited resonant 87 Rb2 molecules by 
collisions with other Rb atoms, based on the results of 
recent experiments |l4j . 

The theory presented in Refs. IE IS El- based 
on coupled Gross-Pitaevskii equations for atomic and 
molecular condensates, cannot be applied to the analy- 
sis of these experiments involving a non-condensed ther- 
mal gas. The approach used here is based on the Breit- 
Wigner theory of resonant multichannel collisions (see 
e.g. Ref. 0), as has been proposed for the system un- 
der consideration by Ref. ■ The reaction involving the 
excited resonant molecule Rb2 (to) includes a reversible 
input channel of formation from (and dissociation to) a 
pair of colliding atoms, 



Rb + Rb 



Rb2 (m) 



(1) 



and irreversible output channels of exoergic collisions 
with a third atom, 



Rb 2 (to) + Rb -► Rb 2 (d) + Rb, 



(2) 



bringing the molecule down to one of the lower-lying rovi- 
brational levels of the same spin state, or to levels be- 
longing to other spin states. (An alternative approach, 



presented in Refs. 0, , treats the whole process as a 
one-stage recombination by a three-body collision.) 

Let us consider all atoms, for the time being, as dis- 
tinguisha ble p articles. According to the standard theory 
(see Ref. [15|), the natural resonance width T e associ- 
ated with channel Q is two times smaller than the cor- 
responding width for the case of indistinguishable atoms 
presented in Ref. (see also Refs. 00). It exhibits a 
Wigner threshold dependence of the form 



r K = 



\a a ii\A 



(3) 



where a a is the non-resonant (background) elastic scat- 
tering length, /i is the difference of the magnetic momenta 
of the atomic pair and the Rb 2 (to) molecule, A is the 
phenomenological resonance strength (see Refs. 00), 
and p is the relative momentum of the colliding atoms. 
These parameters also describe the variation of the elastic 
scattering length function of the external mag- 

netic field B in the vicinity of the resonance at B = Bq 
as (see Refs. 00) 



(4) 



B-B 



The total width Td associated with the deactivation chan- 
nel (J2J can be expressed in terms of a two-body rate co- 
efficient kd, as 



T d = kdn 



(5) 



is proportional to the atomic density n. The rate co- 
efficient kd includes the contributions of all the output 
deactivation channels (d) of Eq. J3J). 

The Breit-Wigner theory leads to the following expres- 
sion for the cross section of resonance-enhanced three- 
body recombination (see Ref. 0), 



irh 2 



r P T ri 



p 2 fj, 2 (B-B y/h 2 + (r e + r d y/4 



(6) 



This expression does not take into account the indistin- 
guishability of the three participating atoms, in which 
case the cross section should be Cmd = 3!<r (see Ref. 15] ). 

The resonant molecular state Rb 2 (m) can be formed 
whenever the detuning from the resonance is compara- 
ble or less than r e . This state decays producing atoms 



2 



with a kinetic energy spectrum of width hT e . Under the 
conditions of the experiments 01 ( a a ~ 98.96 atomic 
units, /i « 2.8 Bohr magnetons, A w 0.17 G for the 
strongest resonance at 1007.34 G in 87 Rb and a collision 
energy of p 2 /m 2/iK) the width calculated with Eq. 
© is given by HT e /kB ~ 7/iK, where ks is the Boltz- 
mann constant. Therefore this energy is less than the 
trap depth of « 20/iK and a spontaneous dissociation of 
the resonance molecule cannot lead to a significant 
loss of trapped atoms (as opposed to the case of a BEC 
- see Ref. @,H|)- Each deactivation event © leads to 
the simultaneous loss of three atoms. Therefore, the loss 
rate for the atomic density n (r, t) can be written in the 
form, 



n (r, t) 



„2p 
3 — cr i: 
m 



ld n 2 (r,t) = -K 3 n 3 (r,t), (7) 



where 



36irh 2 kd\a a [J,\ A 



p 2 (B - B Q ) 2 + H 2 T 2 J4 



(8) 



is the three-body loss rate coefficient. Here the partial 
inelastic width Td is neglected in the denominator in com- 
parison to r e . Even very close to the resonance, as long 
as \B — Bq\ > 0.1 G, the width T e may as well be ne- 
glected, leading to an expression similar to Eq. (9) of Ref. 
[l2l | for the loss in a BEC. However, the rate coefficient 
given by Eq. (JHJ) is six times larger than the correspond- 
ing rate for a BEC. This difference, due to the effects of 
quantum statistics, has been predicted for non-resonant 
three-body recombination in Ref. [T3 |. and observed in 
experiments [l9j . 

In the case of a BEC the atomic density profile is deter- 
mined by the repulsive interaction between atoms. This 
interaction can be neglected whenever its characteristic 
energy, proportional to the elastic scattering length, is 
small compared to the kinetic energy of atoms, 



4tt 



h 2 a rcs n <C k B T. 



(9) 



For the temperature T = 2/iK used in the experiments 
[Tij this condition is obeyed whenever \B — Bq\ > 0.01 G. 
Therefore we can consider the gas as an ideal one with the 
equilibrium density profile described by the Boltzmann 
distribution in the trap potential. 

The loss rate given by Eq. Q is density dependent. In 
the case of an inhomogencous trapped gas the loss pro- 
cesses modify the equilibrium density profile, leading to 
an atomic drift which tends to compensate for this de- 
formation. The characteristic time for this compensation 
can be estimated as the trap period. In the experiments 
[l4| the magnetic field that brings the system close to res- 
onance has been applied during a time interval of t = 50 
ms. This time substantially exceeds the radial trap pe- 
riod (the radial trap frequency is uj r /2ir — 930 Hz), but 
it is less than the axial trap period (the axial trap fre- 
quency is u) a /2ir = 11 Hz). Therefore we can consider the 
radial density profile as an equilibrium one, described by 



a Boltzmann distribution, and write out the atomic den- 
sity profile as 



v(z,t) 
n ( r ' *) = — uT ex P 



x 2 + y 2 



where 



is a non-equilibrium axial profile and 



v{z,t)— I dxdyn(r,t) 



br = ±J 2kBT 



ui r V m 



(10) 



(11) 



(12) 



is the characteristic radius of the atomic cloud. 

Neglecting effects of axial atom transport, a kinetic 
equation for the axial profile can be written in the form 

i>(z,t) = -K 1D v 3 (z,t), K 1D = K 3 / (3irX) . (13) 

The solution of Eq. I|13f) relates the axial profile at time 
t to the initial one at t = as 



i/(r,0) 



y/l + 2K 1D v 2 (r,0)t 



(14) 



Let us suppose that at t — the atoms have a Boltzmann 
distribution with the temperature T and 



i>(r,0) = ^ exp ( ) , 



1A) 



where 



1 

Wo 



2k B T 



(15) 



(16) 



is the characteristic half-length of the atomic cloud and 
iVo is the initial number of atoms. In this case, the num- 
ber of atoms remaining in the trap can be expressed as 



N(t) 



exp (-C 2 ) 



4^ f<% , 
vW Vl + 2 ^i^o* cx P(- 2 C 2 ) 



(17) 



where £ = z/b a . 

Equation lfT7|l . in combination with Eqs. |J3J, ©, |Q 
and {T5)| , allows us to estimate the value of kd by a fit 
to the number of remaining atoms measured in Ref. 0] 
for iVo = 2.8 x 10 6 . The fit produces the optimal value 
of kd = 0.7 x 10~ 10 cm 3 /s. This value is comparable to 
corresponding estimates for Na resonances (1.6 x 10~ 10 
cm 3 /s in Ref. |;4x 10" 10 cm 3 /s and 10" 11 cm 3 /s in 
Ref. [H following the theory of Ref. 0). The results of 
calculations for several values of kd are presented in F ig. [J 
in comparison with the experimental results of Ref. [14| . 

The authors are most grateful to Dr. Stephan Diirr for 
providing a preprint of Ref. 14] and clarifying details of 
the experiment. 
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FIG. 1: Number of remaining atoms as a function of the 
magnetic field in the vicinity of the 1007 G resonance in 87 Rb 
calculated with Eq. (1171 for three values of the deactivation 
rate coefficient, kd = 7x 10~ cm /s (solid line), 10" 10 cm 3 /s 
(long-dashed line), and 5 x 10 -11 cm 3 /s (short-dashed line). 
The circles represent the experimental results of A. Marte et 
al. 13. 
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